Global Monopole in General Relativity 

Kirill A. Bronnikov^f ? Boris E. Meierovich^ and Evgeny R. PodolyakJ 

I Center for Gravitation and Fundamental Metrology, VNIIMS, 3-1 M. Ulyanovoy St., Moscow 117313, Russia; 
Institute of Gravitation and Gosmology, PFUR, 6 Miklukho-Maklaya St., Moscow 117198, Russia 

X P.L. Kapitza Institute of Physics Problems, 2 Kosygina St., Moscow 117334, Russia 

We consider the gravitational properties of a global monopole on the basis of the simplest Higgs scalar triplet model 
in general relativity. We begin with establishing some common features of hedgehog-type solutions with a regular 
center, independent of the choice of the symmetry-breaking potential. There are six types of qualitative behavior 
of the solutions; we show, in particular, that the metric can contain at most one simple horizon. For the standard 
Mexican hat potential, the previously known properties of the solutions are confirmed and some new results are 
obtained. Thus, we show analytically that solutions with monotonically growing Higgs field and finite energy in the 
static region exist only in the interval 1 < 7 < 3, 7 being the squared energy of spontaneous symmetry breaking 
in Planck units. The cosmological properties of these globally regular solutions apparently favor the idea that the 
standard Big Bang might be replaced with a nonsingular static core and a horizon appearing as a result of some 
symmetry-breaking phase transition on the Planck energy scale. 

In addition to the monotonic solutions, we present and analyze a sequence of families of new solutions with 
oscillating Higgs field. These families are parametrized by n , the number of knots of the Higgs field, and exist for 
7 < 7„ = 6/[(2n + l)(n + 2)] ; all such solutions possess a horizon and a singularity beyond it. 



1. Introduction 

According to the Standard cosmological model ^ , the Universe has been expanding and cooling from a split second 
after the Big Bang to the present moment and remained uniform and isotropic overall in doing so. In the process of its 
evolution, the Universe has experienced a chain of phase transitions with spontaneous symmetry breaking, including 
Grand Unification, electroweak phase transition, formation of neutrons and protons from quarks, recombination, 
and so forth. Regions with spontaneously broken symmetry, which are more than the correlation length apart, 
are statistically independent. At interfaces between these regions, so-called topological defects necessarily arise. 
A systematic exposition of the potential role of topological defects in our Universe is provided by Vilenkin and 
Shellard The particular types of defects: domain walls, strings, monopoles, or textures are determined by the 
topological properties of vacuum @. If the vacuum manifold after the breakdown is not shrinkable to a point, then 
the Polyakov-t'Hooft monopole- type solutions |3[ Q appear in quantum field theory. 

Spontaneous symmetry breaking (SSB) plays a fundamental role in modern attempts to construct particle theories. 
A symmetry in this context is not necessarily associated with space-time transformations. It can be a kind of 
"internal" symmetry as well, such as the Grand Unification symmetry, the electroweak and the isotopic symmetry, or 
even supersymmetry, whose transformations mix bosons and fermions. Topological defects, caused by spontaneous 
breaking of internal symmetries (independent of space-time coordinates), are called global. 

A fundamental property of global symmetry violation is the Goldstone degree of freedom. In the monopole case, 
the term related to the Goldstone boson in the energy-momentum tensor decreases rather slowly away from the 
center. As a result, the total energy of a global monopole grows linearly with distance, in other words, diverges. 
Without gravity such a divergence is a general property of spontaneously broken global symmetries. In his pioneering 
paper |^ Polyakov mentioned two possibilities of avoiding this difficulty. The first one was to combine a monopole 
with a Yang-Mills field. This idea was independently considered by t'Hooft Q. This, among other reasons, gave 
rise to numerous papers on gauge (magnetic) monopoles. The second possibility was to consider a bound monopole- 
antimonopole system, whose total energy would be large (proportional to the distance between the components) but 
finite. 

One more opportunity is to take into account the self-gravity of global monopoles, which can in principle remove 
the above self-energy problem. This is also necessary for potential astrophysical applications. Such a study was first 
performed by Barriola and Vilenkin who found that the gravitational field outside a monopole is characterized by 
a solid angle deficit proportional to the SSB energy scale. Harari and Lousto showed that the gravitational mass of 
a global monopole, calculated using the Tolman integral, is negative 0. Solutions with a horizon for supermassive 



global monopoles were found by Liebling |10|, who also confirmed the estimate of Ref. |14 for the upper value of 
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the symmetry breaking energy compatible with a static configuration. The existence of de Sitter cores inside global 
monopoles and other topological defects gave rise to the idea of "topological inflation" |2[ ^ . 

For global strings in flat space, the energy per unit length (without gravitation) also diverges with growing distance 
from the axis, but only logarithmically. However, in general relativity integration over the cross-section yields a flnite 
result 1^, 0]. The gravitational interaction thus leads to self-localization of a global string. Will a similar effect take 
place for a global monopole? An attempt to answer this question, which appears not to be answered in the existing 
papers, was one of motivations for reconsidering the gravitational properties of a global monopole. 

The previous studies have used the boundary condition according to which the symmetry-breaking potential 
should vanish at spatial inflnity. Our approach is different: we do not even assume the existence of a spatial 
asymptotic but require regularity at the center and try to observe the properties of the whole set of global monopole 
solutions. In doind so, among other quantities, we discuss the behavior of the total scalar fleld energy, which turns 
out to be flnite in static regions of supermassive global monopoles. 

In Sec. 2 we present the complete sets of equations for a static spherically symmetric gravitating global monopole 
in two most convenient coordinate systems, namely, with quasiglobal and harmonic radial coordinates. The general 
properties of static global monopoles are summarized in Sec. 3. In Sec. 4 we analyze analytically and numerically the 
speciflc features of a global monopole in the particular case of the "Mexican hat" potential. Sec. 5 contains a general 
discussion of our results, including their possible cosmological interpretation. 

2. Equations and boundary conditions 
2.1. General problem setting 

We begin with the most general form of a static, spherically symmetric metric, without specifying the radial coordinate 
= u: 

ds^ = g^.^dx^'dx" = e^^°dt^ - e^^'du^ - e^^"dn^. (1) 

Here dfi^ — dO^ -I- sin^ Od^p^ is the linear element on a unit sphere and F^^ Fi, and Fn are functions of u . 
The nonzero components of the Ricci tensor are (the prime denotes d/du) 

< = e-'^nF[,' + F;,{-Fi + 2Fl, + F[,)]; 

i?/ = e-^^'[Fl,' + 2F^ + 2Fl,^ + Fl,^-Fi{2Fl, + Fl,)]- 

Ri ^ Ri^- c-2^" + e-2^^ [F^ + Fl,{-F[ + 2Fi + F^)]. (2) 
Consider the Lagrangian describing a triplet of real scalar fields (/)° (a = 1,2,3) in general relativity: 

L = ^ + \g'"'d^<l,'^d,r-V{<i,), (3) 

where R is the scalar curvature, ¥{(()) is a potential depending on cj) — :k\f(\f^ and G is the gravitational constant. 
We use the natural units such that 

;i = c = 1, (4) 

so that G = Tipi^, where mpi = 1.22 x 10^^ GeV is the Planck mass. 

To obtain a global monopole with unit topological charge , let us assume that the metric has the form (|I]) while 
</)° comprise the following "hedgehog" configuration: 

(p"^ — (p{u)n°' , n°' = {sin 6 cos (fi, sin 6* sin cp, cos 6'}. (5) 

The Einstein equations can be written in the form 

r; = -8^Gf; = -8^g(t; - ^s;t^) (6) 

where is the energy-momentum tensor and the nonzero components of are 

fo° = -V, fl = -V ~ e-2-Fi0'2^ ff = ff = -y - e-2f^'^^2 (-y-) 

The regular center conditions for the metric (0) are that at the corresponding value Uc of the coordinate x^ = u 
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The last condition is necessary for local flatness and provides the correct circumference to radius ratio for coordinate 
circles at small r = e^" . 

The scalar field energy, defined as the partial time derivative of the scalar field action, E = —dS/dt, is a conserved 
quantity for our static system: 

"1 



du, (9) 



where g is the metric tensor determinant. 

In what follows we will make some general inferences without specifying the potential V{(j)) and then perform a 
more detailed study for the simplest and most frequently used symmetry-breaking potential 

v{<i>) = ^x{rr - v'? = Iv'Kf 1)', (10) 

where 77 > characterizes the energy of symmetry breaking, A is a dimensionless constant and f{u) = (j){u)/ri is the 
normalized field magnitude playing the role of an order parameter. The model has a global 50(3) symmetry, which 
can be spontaneously broken to SO{2) due to the potential wells (y = 0) at f — ±1. 

Let us now write down the Einstein equations and the boundary conditions explicitly in two coordinate frames 
to be used. 

2.2. The quasiglobal coordinate p 

The first choice is the coordinate u = p specified by the condition Fq + Fi = . Denoting e^^" = e^'^^^ = A{p) and 
g-Fn _ 7^(p)^ -vve obtain the metric in the form 

ds'^A{p)dt'-^^^r\p)dn\ (11) 
The scalar field equation 

act)" + dV/d(t)'' =0, (12) 
where □ = V"Vq is the d'Alembert operator, and certain combinations of the Einstein equations have the form 

(Ar^cj)')' -2(j) ^ r^dV/d(j); (13) 

{A'r'^y = -167rGrV; (14) 

2r"/r = -SttGc/)'^; (15) 

A{r^y' -r^A" = 2{l-87rG(l)'^); (16) 

A'rr + Ar'^ - 1 ^ %'KG[\Ar^(j)''^ - 0^ _ r^y], (17) 

where the prime denotes d/dp. Only three of these five equations are independent: the scalar field equation ( [l3| ) 
follows from the Einstein equations, while Eq. (|l^) is a first integral of the others. Given a potential V{(p) , this is a 
determined set of equations for the unknowns r. A, <f>. 

This choice of the coordinates is preferable for considering Killing horizons, which correspond to zeros of the 



function A{p) , since such zeros are regular points of Eqs. ([13|)-(17); moreover, in a close neighborhood of a horizon, 
the coordinate p defined in this manner varies (up to a positive constant factor) like manifestly well-behaved Kruskal- 
like coordinates used for analytic continuation of the metric JT6| , p^ . Therefore one can jointly consider regions at 
both sides of a horizon in terms of p, and, in general, the whole range of p can contain several horizons. For this 
reason the coordinate p can be called quasiglobal. 

The regular center conditions (||) are fulfilled if, near some value pc of the coordinate p , 

A{p)^A, + 0{{p-p,f), r{p)^{p-p,)/^,. (18) 

In regions where A < Q (sometimes called T-regions), if any, the coordinate p is timelike and t is space-like. 
Changing the notations, t — > a; G K, and introducing the proper time of a comoving observer in a T-region, 



r = j dpl^\A{p)l (19) 

we can rewrite the metric in the form 

ds^ =dT^ -\A{T)\dx^ -r'^{T)d9:^. (20) 

The space-time geometry then corresponds to a homogeneous anisotropic cosmological model of Kantowski-Sachs 
(KS) type where spatial sections have the topology R x §^ . 
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2.3. The harmonic coordinate u 

Another convenient variable, well simplifying the form of the equations, is the harmonic coordinate u specified by 
the condition 

Fi^2Fn + FQ, (21) 

so that □ u = . Then the field equations may be written as 

<j)" -2e^°+^'<j) = e'^^'dV/d(t), (22) 

F;^' = -SnGe^^'V, (23) 

F[' - 2F^{Fl^ + 2F^) = -87rG(<^'2 + e^^^V), (24) 

F/; - e2(^°+^'^) = -87rG(</)2e2(^°+^") + e^^^F), (25) 

- e-2^'^ + e-^^H^/i' + Si^n^o) = 8^G(i e-^^'cj,'^ - e-^^^cf'' - V) (26) 

where the prime denotes d/du. 

It is straightforward to obtain that a regular center can only correspond to u — s- ±oo; we choose u —oo, where 
one should have 

e-^"^ ~ l/\ul = ^/A'e{l + 0{u-^)), e^' ~ (27) 

and Ac is the same as in (p^). 

3. General properties of global monopoles 

3.1. Monopoles in Minkowski space-time 

The Minkowski metric in the usual spherical coordinates 

ds^ = dt^ - dr^ - r^dVi^ (28) 

is a special case of ( [ll| ) with r = p and A = 1. The only unknown in flat space-time is 0(r) , and the only field 
equation is ([l3| ) which takes the form 

(r^^')' ~2(t)^ r^dV/d(t>, (29) 

where, in particular, for the potential ( p^ dV/d(p = X<pi(p — if') - In this case the scalar field equation can be written 
in terms oi f — (jy/rj as 

r-^{r^n' - ^fr-^ + Wfi^ - f) = 0- (30) 
The energy integral (^ takes the form 
^1 



E = Att I 



2 



dr, (31) 



and for its convergence, in case V{(j)) > 0, all the three terms should vanish quickly enough at infinity: 

= o(r-i/2), 0' = o(r-3/2)^ V = o{r-^) (32) 

as r oo . This actually means that a configuration with finite energy is only possible with ^(0) = 0, contrary to 
the symmetry breaking assumption according to which V has minima in nonsymmetric states, <j> 0. In particular, 
the potential (^0|) does not give rise to global monopoles with finite energy. A consideration of self-gravity of the 
field triad is one of the ways to overcome this difficulty. 

The harmonic coordinate u in flat space-time is connected with r by the relation u — uq — ±l/r, where uq is 
an arbitrary constant; choosing the minus sign, we find that u ranges from ~oo, which corresponds to the center 
r = , to uo corresponding to spatial infinity. 

''A cylindrical version of the harmonic radial coordinate has been used previously in the analysis of gravitational properties of current- 
conducting filaments jl^ ] and cosmic strings |l9| pc| . 
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3.2. Solutions with constant cj) 

Under the assumption (j) — (j)o ~ const, the corresponding value of the potential ^(^o) — Vo (times SttG) plays the 
role of a cosmological constant, and the Einstein equations can be integrated explicitly. 

Indeed, in a region where (j) = const, Eq. ( |l5|) reduces to r" ~ 0, whence r = ap + tq , a, rg = const. It remains 
to find A{r) , and this is immediately done by integrating Eq. (^|): 

A{r)^ - + Cr^, A = SttG^^ ^ (33) 

r 

where M and C are integration constants. Substituting ( p3| ) into (p^), we find 

C = -SttGVq/S. (34) 

Thus the solution is essentially determined by the values oi (f>o, Vq and M. One more constant, a, reflects the 
freedom in choosing the unit of time. It should be noted that this is not a monopole solution. Even if we put M = 0, 
which is evidently necessary for regularity at r = 0, this solution with constant 7^ is singular at the center: for 
(H) with M ^ the Kretschmann scalar at small r is /C = Rap-ysR"''^^^ ~ AA'^/r'^. 

With respect to a global monopole, two cases of the solution ( |3^ ) are of interest. The case = describes the 
symmetric state, while the case Vq = gives a possible asymptotic behavior at spatial or temporal infinity. 

In case (po — (the symmetric state), putting M = (which is necessary for a regular center), we arrive at the 
de Sitter metric 



This metric has a horizon at r = ■ At r > , outside the horizon, r becomes a timelike coordinate, and i is a 



spacelike one. Changing the notations as in (19), (|20|), we obtain the metric 



ds^ = - sinh^ {T/rh)dx^ ~ r\ cosh^ {T/rh)dQ'^. (36) 

This is a KS cosmology with isotropic inflationary expansion at late times ( r —^ 00). 

In the other case, (j)o ^ but Vq = (the case of broken symmetry, such as (/) = 77 in the potential (p^), the 
metric takes the form |^ 

ds^ = 5 ]dt^ - { 5 dr^ - r^dn^, (37) 

\ a ^ / \ a T ) 

where the constant M has the meaning of mass in the sense that test particles at rest at large r experience an 
acceleration equal to —GMjr^ in the gravitational field (p7|). Furthermore, a nonzero value of 0o leads to a solid 
angle deficit A defined in ( p3| ) at the asymptotic r 00 (see for more detail) and to a linear divergence of the 
integral (||) at large r. 

The general case of Eq. ( p3[ ) describes the large r asymptotic of any solution to Eqs. (p^-(p!7|) , provided that such 
an asymptotic exists and tends quickly enough to a constant value. 

For the monopoles to be studied, the metric ( p7| ) gives a large r asymptotic in the case A < 1 . We will also 
consider solutions with A > 1, for which a static asymptotic is absent. Then the metric (|37| ) describes cosmological 
evolution at late times. 

3.3. General properties of solutions v^rith varying 

Consider now the general form of Eqs . (|l^ )~(P^) with varying i/i, without specifying the potential V^((/>) . 

First of all, we note that, due to (|15|), r" < 0, which rules out any nonsingular configurations without a center 



such as wormholes and horns (see Theorem 1 of Ref. |17 for more detail). 
Second, Eq. ( [l6| ) can be rewritten in the form 

(r^B')' = -2(1-8^G02)^ Ajr', (38) 

and at a point where 5' = we have r^B" = —2(1 — SttGc/)'^) . Hence it follows that, as long as cfP < 1/{8ttG) (the 
(f) field does not reach trans-Planckian values), B" < at possible extrema of the function B. In other words, B 
cannot have a regular minimum. 

Our interest is in systems with a regular center satisfying the conditions (p^, so that both A{p) > and B{p) > 
near p = pc- At a possible horizon p ~ h, both A and B vanish, and, since it cannot be a minimum of _B, _B < at 
p > h near the horizon. At greater p the function B{p) , having no minima, can only decrease and will never return 
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Figure 1: Possible behavior of r{p) in global monopole solutions 



to zero; hence A = Br^ < at p > h. We conclude that there can be no more than one horizon, and, if it exists, it 
is simple (corresponds to a simple zero of A{p)). Since the global causal structure of space-time is determined (up 
to possible identifications of isometric hypersurfaces) by the number and disposition of Killing horizons j2^, |2^, , 
we have the following result: 

Statement 1. Under the assumption that (f)^ < l/{8nG) in the whole space, our system with a regular center can 
have either no horizon, or one simple horizon, and in the latter case its global structure is the same as that of de 
Sitter space-time. 

The above reasoning is in essence the same as in the proof of Theorem 2 of Ref. on the disposition of horizons 
in scalar-vacuum space-times. It uses only Eq. ( |l6| ) where the potential V does not enter. Therefore the conclusion 
is valid for systems with any potentials, positive or negative. 



Let us return to Eq. (15), according to which r" < 0. Since at a regular center r' > 0, this leaves three possibilities 



for the function r{p) (see Fig. 1): 

(a) monotonic growth with a decreasing slope, but r — > oo as p ^ oo , 

(b) monotonic growth with r — s- r,„ax < oo as p — s- oo , and 

(c) growth up to Tniax at some pi < oo and further decrease, reaching ?' = at some finite P2 > Pi ■ 

In each case, according to Statement 1, a horizon can occur at some p = h within the range of p, so that at 
p > ft, we have a T-region with the geometry of a KS cosmological model. 

We conclude that there are six classes of qualitative behaviors of the solutions, i.e., (a), (b), (c), each with or 
without a horizon, the latter circumstance to be labelled with the symbol 1 or 0, respectively. Thus, all solutions 
with a spatial asymptotic belong to class (aO). Class (bO) includes space-times ending with a "tube" consisting of 
two-dimensional spheres of equal radius. Class (cO) solutions contain a second center at p — P2, and this center 
can a priori be regular or singular. We thus obtain a static analogue of closed cosmologies. Classes (al), (bl), (cl) 
describe different late-time cosmological behaviors in the two directions corresponding to , whereas the fate of the 
third spatial direction (K) is determined by the function A{p) . In particular, the possible de Sitter asymptotic ( |36| ) 
belongs to class (al) solutions, and in this case the expansion is isotropic at late times. On the other hand, class 
(cl) contains models which behave at late times like the Schwarzschild space-time inside the horizon, contracting to 
r = 0. 

This classification is obtained without any assumptions about V{(t>) . Solutions with given V{4i) will contain some 
of these classes, not necessarily all of them. 

In case > , Eq. ( p^ leads to one more important observation: since at a regular center A'r^ = , we can 
write (^4|) in the integral form 

A'r^ = -WttG r V{p) r^{p) dp. (39) 

Thus, unless V = 0, A{p) is a decreasing function. Eq. ( p9| ) leads to the following conclusions. 

Statement la. If V{(f)) > 0, our system with a regular center can have either no horizon, or one simple horizon, 
and in the latter case its global structure is the same as that of de Sitter space-time. 
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Figure 2: The field magnitude / as a function of tlie liarmonic coordinate u for different values of 7. Solutions with 
monotonically growing / = fo{u) (a) exist for < 7 < 70 = 3. In the region 7 < 71 = 2/3 there are solutions with 
/ = fi{u), changing their sign once (b), in the region 7 < 72 = 0.3 there are solutions with / = /2(w) changing 
their sign twice (c). As 7 ^ 7n — 0, the function /„(u) vanishes in the whole range — 00 < u < 00 from the center 
to the horizon. 



Statement 2. If V{<f>) > 0, the second center in class (cO) solutions is singular. 

Statement 3. If V{(f>) > and the solution is asymptotically flat, the mass M of the global monopole is negative. 

Statement la shows that, for nonnegative potentials, the assumption 0^ < l/(87rG) in Statement 1 is unnecessary, 
and the causal structure types are known for any magnitudes of (p. 

Statement 2 follows from A'{p2) < 0, whereas at a regular center it should he A' = 0, see (|l8|). A'{p2) = could 
only be possible with V = 0, but in this case the only solution with a regular center is trivial (flat space, (f> = 0). 

In Statement 3, asymptotic flatness is understood up to the solid angle deficit, i.e., r = p and A is given by ( |33| ) 
with C = at large p. Then, (^9|) for p ^ 00 gives 2GM on the left-hand side and a negative quantity on the right. 

To our knowledge, this simple conclusion, valid for all nonnegative potentials, has been so far obtained only 
numerically for the particular potential (p^ 1^. Note that Statement 3 is an extension to global monopoles of the 
so-called generalized Rosen theorem |2^, |Tj, previously known for scalar- vacuum configurations. 

Thus, even before studying particular solutions with the potential (|l0|), we have a more or less complete knowledge 
of what can be expected from such global monopole systems. 

4. Mexican hat potential 

4.1. Equations and boundary conditions 

Further analysis is performed for the particular "Mexican hat" potential (|l^). For numerical integration we prefer 
to use the harmonic coordinate u and to work with Eqs. (|2^)-(|25|). This variable enters into the equations only via 
derivatives and is thus invariant under translations u u + const . 
Introducing the dimensionless quantities 

u = u/{^/\r|), e-^" = \/Arye^", e^^ = \rj^ e^^ (40) 

we exclude the parameter A from the equations. Indeed, omitting the tildes, we obtain 

/" = c2(^°+^")[2- e2^"(l-/2)]/, (41) 

K = -Je2(^«+2Fn)(j2_^)2^ (42) 

Fli - e2(^»+^n) [1 - 7/2 - J e^^" (1 - ff] (43) 



The condition (|21|) is preserved for the newly defined quantities, but the metric now reads 

2Fi 7 2 I 2Fn^o2 

ds^ = e2^« dt^ - ^ "'^ + ° "'^ . (44) 

The boundary conditions at u ^ — oo are 
/ = 0, Fo = 0, F^ = 0, Fn = - In(-u) + o(l/|u|), as u ^ -oo, (45) 
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They follow from the requirement of regularity at the center and a particular choice of the time unit {Fq = 0) and 
the origin of the u coordinate (the fourth condition). 

There remains only one dimensionless parameter in Eqs. (^)-(^), 

j^SnGr]^, (46) 

it is the squared energy of symmetry breaking in Planck units. 

It is easy to obtain that 7=1 is a critical value of this parameter. Indeed, if we suppose the existence of a large 
r asymptotic at which / ^ 1, i.e., the field tends to the minimum of the potential (|l0[), then the asymptotic form 
of the metric at large r is (|3^) with A = 7 . Consequently, the asymptotic can be static only if 7 < 1 , whereas for 
7 > 1 the large r asymptotic can be only cosmological (KS type), and there is a horizon separating such an outer 
region from the static monopole core. 

On the other hand, if a configuration with 7 < 1 possesses a horizon, there is again a KS cosmology outside it, 
but there cannot be a large r asymptotic and, according to Sec. 3, the solutions belong to classes (bl) or (cl). 

Now, leaving aside the rather well studied case of solutions with a static asymptotic ||, ^ , belonging to class 
(aO) according to Sec. 3, let us suppose that there is a horizon and return to Eqs. (p])-(|43|). The horizon corresponds 
to u — > +00. For such cases, in addition to (E5|), we impose the boundary condition 



/(w)-^A, |A|<oo as -it^+cx). (47) 

This condition is necessary for regularity of a solution on the horizon and is applicable to classes (al), (bl), (cl). 

For class (aO) solutions, having a spatial asymptotic and no horizon, the condition ( p7| ) is meaningless. Moreover, 
the coordinate u then ranges from —00 to some wq < 00 such that r(uo) = 00. 

For configurations of classes (aO) and (al), the commonly used boundary condition at large r is 

/ ^ 1 as r ^ c». (48) 

It is of interest that in the case (al], to which both condition are applicable, the condition (^7|), being less restrictive, 
still leads to solutions satisfying ([48|) due to the properties of the physical system itself. 

comprise a well-posed nonlinear 



The set of equations (^)-(^3|) with the boundary conditions (^) and ( |47 



eigenvalue problem. Its trivial solution, with f — and the de Sitter metric (35|), describes the symmetric state 
(with unbroken symmetry). Nontrivial solutions, describing hedgehog configurations with spontaneously broken 
symmetry, can be found numerically and yield a sequence of eigenvalues 7„ , n — 0, 1, ... and the corresponding 
values of the horizon radius „ for each given value of fh . Conversely, for given (admissible) value of 7 one obtains 
a sequence of values of fh and . 



4.2. Linear eigenvalue problem 

Liebling has found empirically the upper critical value 70 « 3 for the existence of static solutions jl^^. In this 
section we find a theoretical ground for this limit. Actually, we find analytically a sequence of critical values 7„, 
n = 0, 1, such that for 7 < 7n there exist static configurations with the field magnitude f{u) changing its sign n 
times. 

In the literature one can find only an analysis for f{u) > 0. Our numerical integration of Eqs. (|4l|)-(|4^) shows 
that, in addition to solutions with monotonically growing f{u) (which exist for 7 < 70 = 3, see Fig. 2a), for 
7 < 71 = 2/3 there are also regular solutions with f{u) changing sign once, see Fig. 2b. For 7 < 72 = 0.3 there are 
solutions with two zeros of f{u) (Fig. 2c), etc. All these solutions have a horizon, and the absolute value of / on the 
horizon \fh,n\ = |/n(oo)| is a decreasing function of 7, vanishing as 7 — *■ 7„ — 0, see Fig. 3. 

As 7 ^ 7n , the function f{u) vanishes in the whole range of u, and it is this circumstance that allows us to find 
the critical values 7„ analytically. In a close neighborhood of 7„ the field f{u) within the horizon is small, ^ 1, 
so that Eq. (^) r educes to a linear equation with given background functions Fq and Fn , corresponding to the de 
Sitter metric~(p5|) . In terms of the dimensionless spherical radius r , Eq. (^ij) takes the form 



dr 



rH 1 



dr 



(2-r2)/ = 0, (49) 



where = \J\2/^ is the value of r on the horizon. The boundary conditions are 



/ 



= 0, |/(r;,)|<oo. (50) 

r— 



^In the notations of Ref. ^ v^3/(87r). 
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Figure 3: 7 dependence of the values of f{u) on the horizon, „ — 



Nontrivial solutions of (^9|) with these boundary conditions exist for a sequence of eigenvalues 7 = 7n , n 
0, 1, 2, and the corresponding eigenfunctions /n(r) , regular in the interval < r < , are simple polynomials: 



fn{r) ^^ak 

k=0 



2fe+l 



Substituting ( |5l| ) into (49), we find the eigenvalues 

<„ = 2(2n+l)(ji + 2), 7„= ^ 

and the recurrent relation 

(2fc-l)(2fc + 2)-r2 



(n+ l/2)(n + 2)' 



(2A: + l)(2fc + 2) 



/c = 1,2 



, -"7 ■ ■ ■ , 



(51) 



(52) 



(53) 



allowing one to express all Ofc, k — 1,2, ...,n, in terms of oq. Eg. (^) is linear and homogeneous, so ao is an 
arbitrary constant^. For fixed n the coefficients in ( |5l| ) are 



^ (2z-l)(2z + 2)-r^„ 
gfc = ao .w„. n > 0, 



(22 + l)(22 + 2) 



- 2 



1< k <n. 



(54) 



The case 
n = 0, 



rh,o 



fo{r) = aor/rtfl 



gives a monotonically growing function /(u) in a close vicinity of 7 = 70 = 3, see Fig. 2a. Thus the upper limit 
7o = 3 for the existence of static monopolc solutions, previously found numerically by Liebling ]lO{ , is now obtained 
analytically. 
The case 



n = 1, 



1 = 3V2, 



fi{r) = ao- 



r 






ao 











describes the function f{u) , changing its sign once, at 7 close to 71 = 2/3, see Fig. 2b. The case n — 2, 72 = 3/10, 
r,,2 = 2^10, 



f2{r) = ao 

rh,2 



5 \rh,2 



99/ r 
35 \rh,2 



gives the field function f(u) changing its sign twice (Fig. 2c). 



^At 7 — > 7n — 0, the general equation ( |4l[ ) has the same solution as (^9|), with ao 1 ■ To find the dependence ao{-f) one has to take 
into account the next terms nonlinear in / . 



10 



/ 



0.0010 



0.0005 



0.0000 



-0.0005 



y= 0.001 




Figure 4: The field magnitude / as a function of the spherical radius r for 7 = 0.001. 



For n » 1 the function /„(r) rapidly oscillates: 



cos 



fn{r) = ao- 



rh^n arcsm , 



^ - 2 . 
-^^ V 2 arcsm 




1 - {r/rh ,n) 
1 - 2/r 



2 \ 1 



4/^2(^2 _ 2) [l-(r/r„,„)2 



However, this semiclassical formula is not valid near the left turning point^ r = \/2 , see dashed curve on Fig. 4. Its 
applicability range is 1 ^ r < r^^^ « 2n , n ^ 1 . 

We have not met so far regular monopole configurations with the field function f{u) changing its sign. It seems 
that this is their first presentation. 



4.3. Solutions with monotonically growing f{u) 

As is clear from the aforesaid, the interval < 7 < 3 of existence of nontrivial solutions with monotonically growing 
f{u) splits into two qualitatively difi^erent regions, separated by 7 = 1. 

In the interval < 7 < 1 the solutions have the spatial asymptotic (^7|) and, according to our general classification, 
belong to class (aO). The spherical radius r{u) = e^"*-"-* varies from zero to infinity, /(u) grows from zero to unity, 
A{u) decreases from unity to its limiting positive value (compare with (37)) 



dr 
dp 



(55) 



and the energy integral (^ diverges. 

In the interval 1 < 7 < 3, the solutions with monotonically growing f{u) belong to class (al). Instead of a 
spatial asymptotic, there is a horizon and a KS cosmology outside it. The functions A and r inside and outside the 
horizon are presented in Fig. 5 for 7 = 2. 

In the presence of a global monopole the cosmological expansion is slower than the de Sitter one (^) . As r ^ 00 , 
the radius r(r) grows linearly, while A tends to the negative constant value —(7 — l)/a^. 

Within the horizon f{u) monotonically grows from zero at m = —00 to a value fh = fh,o{l) on the horizon, 
u — !■ 00 , see Fig. 2a. The / value on the horizon fh.o as a function of 7 decreases from unity at 7 = 1 to zero at 
7 = 7o = 3 , see Fig. 3. The integral (||) taken over the static region converges, and we can conclude that at 1 < 7 < 3 
the gravitational field is strong enough to suppress the Goldstone divergence and to localize the monopole. At 7 > 3 
gravity is probably so strong that it restores the high symmetry of the system. 

Outside the horizon the field / as a function of the proper time r grows from //j q on the horizon to unity at 
T ^ CO. Introducing the proper radial length / inside the horizon by the relation dl — dp/^/A, one can ascertain 
that the functions f{l{p)) at p < h and f{T{p)) dX p > h are two parts of a single smooth curve, see Fig. 6. 

^Recall that in view of the substitution (MOl) the distances are measured in the units {■\/~\ri)~^ . 



11 




Figure 5: The functions A{p) and r{p) form unified smooth curves in the regions inside (sohd curves) and outside 
(dotted) the horizon; 7 = 2. 

/(P) : 




2 p 4 6 8 10 P 



Figure 6: f{p) inside (sohd) and outside (dotted) the horizon; 7 = 2. 



When the parameter 7 is close to its critical value (7 = 1), separating the (aO) and (al) branches of the solution, 
i.e., for 

0<7-l<l. (56) 

one can find analytically the horizon radius r^fi and the scalar field value on the horizon ff^ Q under certain additional 
assumptions on the system behavior which follow from the results of numerical analysis. In particular, there is an 
"intermediate" region of the u range, 1 -C u <C uo = const, where the first term /" in the scalar field equation ( p| ) 
is very small whereas the function e^(-'^''+^") is quite large (despite the fact that this function eventually vanishes as 
u — > 00). Therefore in this region the expression in square brackets should be small, i.e., 

e2^-(l-/')«2, 

and this relation can be used for further estimates. 
The results are 

lnr,,.,o«ln[l/(7-l)]» 1, A,o « 1 - ^(7 - 1)^, (57) 

where the constant C can be found by comparison with the numerical results; our estimate is C w 0.3 

The solution behavior in the critical regime, 7 = 1, can be characterized as a globally static model with a "horizon 

at infinity" since A ^ as r — > 00 . 

The fact that monotonic solutions with horizons are absent for 7 < 1 becomes clear from an analysis of the 

inflection point u = Ujnf of the function Fq{u) . A horizon, if any, corresponds to m — s- cx) where Fq remains finite. 
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-2 
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2.99 



Figure 7: The function F^iiu) for different values of 7. At 7 < 1 there is a Umiting value Umax of u, such that 
F[2 ^ 00 as u itmax(7) (dotted curves). As 7 — > 1, the value Umax(7) 00 (dashed curve), and, for 7 > 1, 
Fn{u) tends to a finite constant value as u — ^- cxd (solid curves). 

Since it behaves logarithmically as u — > —00 , there is (at least one) inflection point, where the second-order derivative 
is zero, and from Eq. (Bs) we have 



A monotonically growing function f{u) corresponds to greater values of /, i.e., to the "minus" branch of ( |5q ) (as is 
confirmed by numerical results). But then the r.h.s. of (^ ) is negative for 7 < 1, leading to > 1, which cannot 
happen since = 1 is the maximum attainable value for the solutions under study. Thus, for 7 < 1 , all solutions 
with monotonically growing f{u) belong to class (aO), possess a spatial asymptotic with a solid angle deficit and a 
divergent field energy. 

Numerical integration confirms these conclusions. The different behavior of Fq{u) for 7 < 1 and 7 > 1 is shown 
in Fig. 7. 

4.4. Solutions with f{u) changing its sign 

For 7 < 71 = 2/3 there are solutions with the function /(0) changing its sign once, see Fig. 2b. For 7 < 72 = 0.3 
there are solutions where f{(j>) changes its sign twice, see Fig. 2c, etc. Unlike the monotonic solutions discussed 
in Sec. 4.3, all of them possess a horizon, and, in agreement with the general inferences of Sec. 4.1, they belong to 
class (cl). This means that, beginning with a regular center, the spherical radius r{p) first grows, then passes its 
maximum rmax at some pi and then decreases to zero at finite p = P2 which is a singularity. The horizon occurs at 
some p = h < p2, which can be greater or smaller than pi , but in any case the singularity takes place in a T-region 
and is of cosmological nature. The dependence r{p) before and after the horizon is a single smooth curve, see Fig. 8a. 

Beyond the horizon, |A(t)| as a function of the proper time r of a comoving observer grows from zero at r = 
(the horizon) to infinity at r ^ = r(p2) (the singularity), see Fig. 8b. The scalar field magnitude |/(r)| beyond 
the horizon first grows, then slightly varies around unity. Approaching the singularity, /(r) changes its sign and 
finally |/(t)| ^00 as r ^ , see Fig. 8c. 

5. Conclusion and discussion 

We have performed a general study of the properties of static global monopoles in general relativity. We have shown 
that, independently of the shape of the symmetry breaking potential, the metric can contain either no horizon, or 
one simple horizon, and in the latter case the space-time global structure is the same as that of de Sitter space-time. 
Outside the horizon the geometry corresponds to homogeneous anisotropic cosmological models of KS type, where 




This is a quadratic equation with respect to 1 — / 



^ , whence 




(58) 
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Figure 8: Solutions with / of alternating sign: the functions (a) r{p), (b) |^(p)|, and (c) f{p) for 7 — 0.5 before 
(solid) and after (dotted) the horizon. 



spatial sections have the topology R x §^ . In general, all possible solutions can be divided into six classes with 
different qualitative behavior. This classification is obtained without any assumptions about T^(0) . Solutions with 
given V{(j)) contain some of these classes, not necessarily all of them. This qualitative analysis gives a complete 
picture of what can be expected for global monopole systems with particular symmetry breaking potentials. 

Our analytical and numerical analysis for the particular case of "Mexican hat" potential confirmed the previous 
results of other authors concerning the configurations with monotonically growing order parameter. Among other 
things, we have obtained analytically the upper limit 70 = 3 for the existence of static monopole solutions, previously 
found numerically by Liebling . We have also found and analyzed a new family of solutions with the field function 
/ changing its sign, which we have not met in the existing literature. 

Of particular interest can be the class (al) solutions with a static nonsingular monopole core and a KS cosmological 
model outside the horizon. Its anisotropic evolution is determined by the functions of the proper time |A(t)| (the 
squared scale factor in the R direction, r^r) (the scale factor in the two directions) and the field magnitude /(r) . 
For a comoving observer in the T-region, the expansion starts with a rapid growth of |^(t)| from zero to finite values, 
resembling inflation, and ending with A — > const as r — > 00 . The expansion in the §^ directions, described by r(T) , 



is comparatively uniform and linear at late times, i.e., much slower than the de Sitter's cosh^(r/r?i) , see (35). It 
should be stressed that all such models with de Sitter- like causal structure, i.e., a static core and expansion beyond 
a horizon, drastically differ from standard Big Bang models in that the expansion starts from a nonsingular surface, 
and cosmological comoving observers can receive information in the form of particles and light quanta from the static 
region, situated in the absolute past with respect to them. Moreover, in our case the static core is nonsingular, and 
it is thus an example of an entirely nonsingular cosmology in the spirit of papers by Gliner and Dymnikova |2^, |29| . 

The nonzero symmetry-breaking potential plays the role of a time-dependent cosmological constant, a kind of 
hidden vacuum matter. Since the field function / tends to unity as t ^ 00, the potential vanishes, and the "hidden 
vacuum matter" disappears. 

The lack of isotropization at late times does not seem to be a fatal shortcoming of the model for two reasons. First, 
if the model is applied for describing the near-Planck epoch of the Universe evolution, then, on the next stage, the 
anisotropy can probably be damped by diverse particle creation, and the further stages with lower energy densities 
may conform to the standard picture (with possible further phase transitions). Second, if we add a comparatively 
small positive quantity A to the potential ( p^ ("slightly raise the Mexican hat"), this must change nothing but the 
late-time asymptotic which will become de Sitter, corresponding to the cosmological constant A. In our view, these 
ideas deserve a further study. 

Evidently, the present simple model cannot be directly applied to our Universe. It would be too naive to expect 
that a macroscopic description based on a simple toy model of a global monopole with only one dimensionless 
parameter 7 can explain the whole variety of early-Universe phenomena. Nevertheless, it may be considered as an 
argument in favor of the idea that the standard Big Bang might be replaced with a nonsingular static core and a 
horizon appearing as a result of some symmetry-breaking phase transition on the Planck energy scale. 
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